Quantum manifestation of Levy— type flights in a chaotic system 
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Semi-classical dynamics of quantum wave packets spreading is studied for a kicked rotor. Quan- 
tum flights are established for a specific, "magic" value of a chaos control parameter when the 
classical stickiness of trajectories is most effective. By studying of survival probability and distribu- 
tion of the accelerations we identify the presence of quantum Levy-type flights. 
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Levy or Levy- type processes occur in different processes related to chaotic dynamics (see reviews in jl]-||]). I n 
Hamiltonian dynamics a phenomenon of stickiness of the chaotic trajectories to island boundaries reveals an intermit- 
tent kinetics which can be partly described by Levy-type processes and fractional kinetic equations . One can find 
more applications in recent publications j^JIIJ . Sticking islands can emerge from a tangent bifurcation [ p~T|JT2) or from 
periodic solutions of parabolic type |]il| , [l4| . The phenomenon of stickiness to islands reveals in power dependence 
of the correlation functions, Poincare recurrences distribution, exit time distribution, etc. 0, O , |la -[l7j . Algebraic 
dependence of these distributions and moments is not universal; it depends on the models of intermittent dynamics, 
time intervals, domain of parameters, etc. Pp^], and even using a single characteristic exponent for the algebraic 
dependence is only a rough approximation to a much more complicated multi-fractal structure of dynamics in chaotic 
sea in the presence of islands. 

Despite the complexity of chaotic dynamics with sticky islands, anomalous diffusion in the form (|p| a ) ~ t" , with 
momentum p and appropriate exponents a and /?, can often be a good approximation, at least under specified 
conditions. The case when fi = 2j3/a > 1 is called super-diffusion, while the case a = 2, j3 = 1 corresponds to a 
normal diffusion. There are different ways to explain the emergence of super-diffusion. One of them, that works for 
the standard map or the web map, is that a trajectory sticks to an accelerator mode island boundary and travels for 
a while almost in the same way as a regular (non-chaotic) trajectory of the accelerator mode |§| |Tl|JT5|Jl7|| . As a result, 
one can observe long flights (Levy-type flights) in phase space, which corresponds to almost regular pieces of chaotic 
trajectories. More precise definition of flights refers to the parts of a trajectory with a Lyapunov exponent A fairly 
small comparing to a typical value of A far from the island boundary layers. 

In one way or another, the phenomenon of Levy-type flights, algebraic correlation decay, stickiness, etc., became 
objects of intensive studies in quantum physics: Levy-type processes of scattering and random matrix distributions 
[ fi"8| , [i"9| recoils distribution in the laser atoms cooling [ pOpl] ] ; anomalous transport in the quantum kicked rotor p3 |2j| 
and in the Harper model [p6| ; algebraic decay of the survival probability in different systems |27-29| and in microwave 
ionization of Rydberg atoms |^,^l| ; fractal conductance fluctuations |H| . 

In this paper we consider a quantum kicked rotor (QKR) system, in which the influence of the sticking dynamics 
on quantum localization process was already demonstrated in |p3 for a specific value of the chaos control parameter 
K = K* = 6.908745 . . ., called "magic" number. It has been shown in |L5| that the phenomenon of stickiness of 
a trajectory to the hierarchical islands-around-islands structure is the most efficient at this magic value of K. The 
trajectories from the chaotic sea can travel over long distances in the phase space together with an island if the 
structure of the boundary islands chain (BIC) is sticky. The stickiness of a trajectory to the hierarchical islands- 
around-islands structure is intimately related to magnitude of the parameter K . It has been shown in p5| that the 
effect of non-Gaussian chaotic dynamics (Levy-type process) becomes more observable when the BIC appears for 
special magic values of K . In particular, for the standard map the value of K = K* = 6.908745 corresponds to the 
sequence of islands 3 — 8 — 8 — 8 — .... This exact self-similarity of the islands structure is also reflected in the scaling 
of the spatio-temporal characteristics of the every island chain generation. A renormalization group was introduced 
in the form T" — > XtT, S' — > XsS (At > 1, As < 1), where S, S' are the areas and T, T' are the periods of the last 
invariant curves inside of the islands, and these characteristics are related to the islands chain and its next generation. 
The scaling constants At and As determine the transport exponent fi = j^j^j ■ , and this theoretical prediction is in a 
good agreement with numerics Jig ]. It should be noted that there are no strict conditions on the hierarchy of islands 
in the case of anomalous transport, and the sequences of values for the islands-around-islands periods can be any 
numbers permitted by the possible resonances of the system. Each magic value of K generates its own sequence of 
islands (i.e., it is 3-8-8-8-. .. for K = 6.908745 .. . |J], or 8-8-8-. .. for K = 6.35933 ... and 5-11-11-11-. . . 
for K = 6.476939 . . . fl?}] ) and its own scaling constants As and At- Therefore different hierarchical sequences specify 
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different deviations from the Gaussian transport, i.e. normal diffusion, and different characteristic time scales when 
the anomalous properties of the diffusion can be observed in both classical and quantum cases. 

An essential increase in the level occupation amplitudes |/ n | of the order of 10 5 has been found for the magic 
value K* for strong quantum (24) and semi-classical [^5| cases. Anomalous diffusion (23) and the quantum tunneling 
from the acceleration mode islands |^| for QKR have been studied as well. These effects are implicit indications of 
Levy-type flights taking place in the quantum system due to the classical sticking dynamics. Accelerator mode islands 
move a quantum wave packet in the energy levels space of the unperturbed system faster than in the case of global 
chaos, since the Lyapunov exponents in the vicinity of these islands are considerably smaller than in the chaotic sea. 
Hence a trajectory sticks to this hierarchical area for long time. The question of direct observation of these flights is 
still open, and we consider this problem here. 

We study numerically the anomalous transport for the Hamiltonian of the QKR 

V 2 

H = — +KcosqS T (t), (1) 

where p = —ih-^ is the angular momentum and q £ [0, 2tt] is the angle, while h is the dimensionless Planck constant, 
K is the control parameter, and Sx(t) = S(t — n) is a periodic train of 8 — kicks with the period T = 1. The 
quantum map 

V(t + T) = U(T)V(t) (2) 

is produced by the following evolution operator over the period T : 

U = cxp{-ihh 2 /2} exp{~i(K/h) cos<?}. (3) 

Below we study the evolution of the complex amplitudes for the expansion of the wave function $ in terms of the 
unperturbed system basis e mq , 

oo 

*(t) = *(g,t)= J2 /»(*)«*"* W 

n— — oc 

Here j/n(i)| 2 is the probability of level occupation at time t which determines the spreading of the energy growth 
over the unperturbed spectrum. It enables us to study the evolution of survival or staying probability in an energy 
region with momentum range Ap — hAN, which is not necessarily compact, and determined without constraints on 
the boundaries: 

n£AN 

For iteration of the map (^) we use the standard technique of the fast Fourier transform with up to 2 18 angular 
momentum eigenstates. We study also the time evolution of the second moment a?. = {p 2 ) = h 2 J2 n \fn\ 2 n 2 as a 
function of K and the initial conditions. The initial conditions for the amplitudes of level occupation are f n (t = 0) = 
<5n,n' where n' = corresponds to the initial population inside acceleration islands, while n' = ±31 correspond to the 
initial population in chaotic sea with p fts i 3 ^ for h — 0.0396. The values of h are irrational numbers or rational 
approximations of irrational numbers such that quantum resonances with 02 ~ t 2 are avoided. 

We started the study of the wave packet spreading for f n (t = 0) = S ni Q. When the chaos parameter K > 1 
approaches the magic number K* — 6.908745 . . . the anomalous super-diffusion takes place in the energy space. In 
this case, sticking dynamics inside of the hierarchical islands around islands structure spreads the initial wave packet 
very fast. This leads to the anomalously diffusive growth of the second moment 02 (i) ~ t@, where (5 f» 1.57 ± 0.27 in 
the time interval t £ [0, 300]. Then the transport exponent changes to diffusive one f3 w 0.98 ± 0.04 determined in the 
time interval t £ [400, 1000] (see Fig. 1). The diffusion coefficient D ~ 10 3 is anomalously large and differs essentially 
from normal diffusion coefficient K 2 /4 <C D. 

The super-diffusion for K = K* is also determined by an increase of |/„| by more than 10 10 times in comparison with 
\f n \ corresponding to usual diffusion with K = 9.69 (see Fig. 2). This differences can be explained by quantum flights, 
which reflect Levy-type flights in the classical limit. To illuminate these flights we compute survival probability in 
periodic energy level regions taken in periodic classical cells of the size 2w x 2tt. To that end ^ = 157 ^ gm is determined 
by a rational approximation of the golden mean value: gm = 55/34, and one obtains for the dimensionless Planck 
constant h = 2ir ■ R/Q = 2ir ■ 34/5393. Therefore h is commensurate with energy level spacing with the period of 
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5393 levels, which corresponds to the 34 cells of the size 2n on a cylinder. One takes 40 levels in every periodic cell 
and studies the change of P(t) in (H) with time, where 

a 

AN = (J [1-Q + sgn(l),l ■ Q + sgn{l) ■ 40] (6) 

l — — a 

and a is an integer part of the ratio N/2Q. The result of the numerical calculation is shown in Fig. 3. Every passage 
of the wave packet through the periodic cell contributes to the increase of survival probability P(t). This causes the 
spikes of P(t). The period between any two spikes is equal to the number of periodic cells: At — R = 34. This means 
that the part of the wave packet which sticks to the BIC passes through one cell per one iteration, i.e. it moves 
ballistically. Amplitudes of the spikes decay exponentially with time by the following expression ~ exp(— 0.005i). It 
relates to the quantum tunneling through cantori [ p2p^j3^ ] . This behavior of quantum tunneling from the accelerator 
mode to the chaotic sea is in a good agreement with the result obtained in |2^ , ^5|] . 

The spikes disappear if the initial conditions are chosen inside of the chaotic sea with n' = ±31 for f n (t = 0) = 8 n . n i. 
Nevertheless, super-diffusion due to the sticking to BIC dynamics takes place for K = K* with 02 ~ t@, where 
j3 ~ 1.39 ± 0.16 . Also, in comparison with K = 9.69, an increase of order of magnitude in the level occupation 
amplitudes \f n \ is obtained. It is the same order as shown in Fig. 2. These differences unambiguously indicate the 
presence of flights. To pinpoint these quantum flights, we calculate a set of the survival probabilities Pk{t) constructed 
over a set of inclosed compact regions of the energy levels AN(k) = Nok, where Nq = 3 6 and k = 1 -j- 10 are integer 
numbers. Let us denote the time intervals which are necessary for the wave packet to reach the boundaries of AN(k) 
as tk- In our numerical calculations tk is the number of iterations. We found a linear growth of the times tk which 
is proportional to the increase of size of the regions AN(k), such that tk ~ AN(k) (see Fig. 4(c)). It means that 
the part of the wave packet which reaches the boundary of AN(k) at the moment tk moves like in an acceleration 
mode: Ap ~ hAN ~ At. As far as the diffusive spreading of the wave packet for K = 9.69 goes, one can see in Fig. 
4(a) that in this case tk is not proportional to AN(k). Another important characteristic of the wave packet evolution 
is the decay rate of the survival probability. Starting from tk the survival probability Pk(t > tk) in (|^) decays in 
correspondence with the rate of the spreading of the wave packet. The decay rate depends on A". In the chaotic 
case when K = 9.69 only the tail of the wave packet moves ballistically. It corresponds to the transitions between 
any two levels n\ and n 2 such that \ri\ — n 2 \ ~ K/h. In this case the matrix elements of the evolution operator 



(||), U ni . ri2 — e m i h / 2 J ni _„ 2 (K/h) ~ yh/K, where J m (z) is the Bessel function of the first kind, are small and 
uncorrelated random quantities. Their total contribution to the transition is summing up to an exponentially small 
value, which is reflected by the exponentially small values of |/„| in Fig. 2(a). Therefore, Pk(t) decay gradually and 
slowly with time as it is shown in 4(a,c). For K = K* the values of Pk(t) fall down abruptly e.i. much faster then 
in the chaotic case, which is illustrated in Fig. 4(b,d). It means that there is a part of the wave packet that moves 
ballistically together with the tail. Obviously that this part determines the decay of Pk(t). Therefore, in case of mixed 
phase space the hierarchical sticking islands structure leads to ballistic flights. Fast change of Pk(t) in Figs. 4(b,d) 
relates to these quantum flights, while the slow changes of Pk(t) in Figs. 4(a,c) relate to normal diffusion. 

It should be noted that the classical nature of the phase space reflects also in pure quantum characteristics of the 
evolution operator U ([|). Let us show that the matrix elements corresponding to flights, such that |ni — r^2 ~ K/h 
are much stronger correlated for K 2tt than for K — 9.69. We consider the following correlation function: 

P(v) = ^2({ U m,nU*,m+r,))m' ( 7 ) 

71 

where {{. . . )) = j; X)m=i ' ' ■ means averaging over m in some energy levels region, such that 1 <C L <C 1/h. By 
using the explicit form of the matrix elements Um^ n , one obtains for the correlation function 

i _ p ihr)L i 

P(V) = je-^' /2 U2K/h)- T -. (8) 

Let us take r\ = [k/h] = K/h — {K/h}, where [Z] is the integer part, while {Z} is the fractional part of Z. Therefore 
hrj = K — h{K/h} = K — e, and the condition Le <C 1 is also fulfilled. Hence, it is easy to see from (@), that if 
K = 2tt + e (or K = 2nn + e), then p(rf) = p a cc{[K/h]) ~ J, K ,^(2[K/h]) is the correlation function with presence 
the acceleration mode motion in the classical counterpart. It is also easy to see for the global chaos case (including 
K = 9.69), that the correlation function satisfies p(rf) = Pdiff([K/h] = i^p acc ([K/h]) <§; p a cc([K/h]). It means that 
for the diffusive case without any stability islands the matrix elements are fairly less correlated than in the case of 
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anomalous diffusion, in which the acceleration mode islands (with K = 2-kji + e) play crucial role for both classical 
and quantum transport and quantum localization of anomalous diffusion |p5| . It should be noted that for K — 2tt 
the motion is diffusive both for classical and quantum case. Only for 2ir < K < 7 the anomalous diffusion has been 
found. Similarly to our result, the anomalous diffusion for QKR has been found in |2|] and [^3| for K = 2n • 1.03 
and K = 6.9115 correspondingly. This range of K relates to bifurcate emerging and existing of the acceleration mode 
islands. That fact that the quantum correlation functions (0) reflects the sticking dynamics in the vicinity of the 
hierarchical structure of these islands relates to the fundamental correspondence principle by which the semi-classical 
wave packet evolves alone an individual trajectory during the quantum correspondence time Th f^J^J^ojl . The h- 

scaling for this time is algebraic, Th ~ h 1 ^ pTj] for the acceleration mode dynamics. This result corresponds to 
general consideration of quantum dynamics with mixed phase space J2^|38| 

To complete this analysis of quantum flights, we consider the acceleration statistics. These values are defined as 
the mean square momentum changes at every iteration, 



where (. . .} means the quantum average. It can be expressed by the level occupation amplitudes in the following form: 



One should note that the dynamical consideration is restricted in time because of the fast spreading of the wave 
packet. This leads to the fast overfill of energy levels region with the fixed size of N = 2 18 levels. Therefore, the 
maximal number of iterations is restricted and equals to 1000 for h — 0.0399. Dynamical calculations of f n (i) are 
carried out for a set of initial conditions, arbitrary chosen inside the chaotic sea to obtain an appropriate statistics. 

Let us consider probability distribution function Wk(Sp) of accelerations Sp. For the white noise Wk(Sp) is 
Gaussian, while for a process which corresponds to the fractional kinetics, one expects Wic(5p) ~ l/(Sp) , where 
\Sp\ — > oo. The exponent I — log 1Q (WK(Sp)/ log 10 (<5p) determines to what kind of statistics Wk(Sp) belongs. When 
I is small enough, we associate the distribution Wx(Sp) with the Levy-type probability distribution function, and 
vice versa, if I is large, Wk (Sp) is associated with the Gaussian one, since one cannot distinguish between large I and 
exponentially small values of Wk (Sp) . The distribution of accelerations Ap = Sp — Sp for the set of 36 different initial 
conditions with 31 < \uq\ < 70, arbitrary chosen inside the chaotic region is plotted in Fig. 5. Here Sp is the mean 
value averaged over t. 

The probability distribution functions are asymmetric and we compare only the right tails of the chaotic dynamics 
with K = 9.69 and the sticky case with K — K* . The distribution function for chaotic, non-sticky case has a short 
right tail with a large slope of I « 9.5. This magnitude of I is large enough, and one concludes that the semi-classical 
diffusion gives rise to Gaussian distribution for the accelerations. The distribution function for K = K* differs 
essentially from the Gaussian distribution. As it is shown in Fig. 5, it is narrower with a longer tail and the slope 
is of the order of I s=s 4.3. These differences enable us to conclude that Wk=k* (Ap) is the Levy- type distribution 
and relates to the quantum flights which are due to the effect of stickiness of the semi-classical wave packet to the 
classical hierarchical islands-around-islands structure. 
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FIGURE CAPTIONS 

Fig. 1 Time evolution of the second moment er 2 for initial population inside the acceleration mode island: n' = 0, 
and h = f^, where (a) K = 9.69, and (b) K = K* . 

Fig. 2 Level occupation amplitudes |/„| vs n for K = 9.69 (a) and K — K* (b) after 1000 iterations. The slop of the 
linear plot for the energy region (6 -j- 12) x 10 4 is 1.4 • 10~ 5 . 

Fig. 3 Semi-log plot of the survival probability P(t) vs t for periodic cells with the period 647r. The plots are: (a) 
for K = 9.69, and (b) for K = K* . 

Fig. 4 A set of survival probabilities Pk{t) vs t for K — 9.69 (a,c) and K = K* (b,d) for variety of energy regions 
AN(k) — 3 6 • k and k = 1 -j- 10 and /i = 2ir/ (157 + .gm) Plots (c) and (d) zoom in of (a) and (b) respectively for 
k = 3. 

Fig. 5 (upper plots) Probability distribution functions Wk(^p) for 120 boxes and the centered by Ap = 5p — 5p, 
where Sp is the mean value with K = 9.69 for (a); and K = K* for (b). (bottom plots) Slops for the right 
tails of the probability distribution functions, where (a) corresponds to K = 9.69 and the slop is 9.51 ± 0.02 for 
Ap e [0.0111, 0.0146]; (b) corresponds to K = K * and the slop is 4.31 ± 0.06 for Ap E [0.0109, 0.0156]. 
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